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Abstract 



It is shown how to adapt the non-perturbative coupled cluster method of many-body 



theory so that it may be successfully applied to Hamiltonian lattice SU(N) gauge theories. 
The procedure involves first writing the wavefunctions for the vacuum and excited states 
in terms of linked clusters of gauge invariant excitations of the strong coupling vacuum. 
The fundamental approximation scheme then consists of i) a truncation of the infinite set 
of clusters in the wavefunctions according to their geometric size, with all larger clusters 
appearing in the Schrodinger equations simply discarded, ii) an expansion of the truncated 
wavefunctions in terms of the remaining clusters rearranged, or "shifted" , to describe gauge 
invariant fluctuations about their vacuum expectation values. The resulting non-linear 
truncated Schrodinger equations are then solved self-consistently and exactly. Results are 
presented for the case of SU (2) in d = 3 space-time dimensions. 



This letter reports the first results of a new approach to Hamiltonian lattice gauge 
theories. The goal is to find approximate non-perturbative solutions of the Schrodinger 
equation for both vacuum and excited states. Space is treated as a discrete lattice in order 
to i) reduce field theory to n — ► oo body quantum mechanics, ii) render the theory ultra- 
violet finite, and iii) make trivial the imposition of Gauss's law, which must be imposed 
as a subsidiary condition on the states in an Aq = gauge, which is the natural choice 
with a Hamiltonian approach]]. 

The basic variables are the SU(N) matrices U{L) defined on each link L. Conjugate 
"chromoelectric" fields E a (L), a = 1 . . . N 2 — 1, are given byQ 



E«(L) - -(l*U(L)) (jW ^- + (v\L)T-) itW ^- j (1) 



in Schrodinger representation, where T a are the group generators. The Hamiltonian for 
SU{N) in d space-time dimensions, without quarks, may be written |2j 

H[x) = a^x'^Wix) (2) 

where the dimensionless variable x is related to the coupling g by x = 4g~ 4 a 2c(_8 , and 

W(x) =Y / E a (L) 2 -xJ2^U p (3) 

L P 

where U p is the trace of the product of Us and U^s around the four sides of a plaquette 
p: U p = TrC/(l)C/(2)C/^(3)C/^(4). The Schrodinger equation for the "reduced" energy 
w = ax l l 2 E is then 

= (4) 

The continuum limit occurs at x — > oo. 

We consider first perturbation theory about the strong coupling x = vacuum, but 
emphasize that the method we shall develop applies at all couplings < x < x c , where 



1 To avoid ghosts, a non-covariant gauge is necessary, which must be Ao = to retain rotational 

invariance. In contrast to QED, Gauss's law cannot be solved explicitly Q. 

2 This operation, in which the elements of the unitary matrices U(L) are treated as if they were 

independent, and with the indices i,j summed over, correctly reproduces the fundamental commutation 
relations for E a (L) with U(L'), U\L') and E a '{L'). 
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x c is the lowest x at which there is a phase transition. For SU(N), it is expected that 
x c = oo i.e. the continuum limit. At x <C 1, the vacuum wavefunction of, for example, 
the SU (2) theory in d = 3 is 



+ * 2 (>En En -^Eo -iirE(m+B) + 

^ s s ' s s s 

+ - 3 fTfeEDEnEn-2 1 6EoEn+----6iioEn) + --- <*> 



where the sums are over all sites s on the lattice. As shown originally by Hubbard S in 
the case of continuum field theory, and observed by Greensite for lattice theories Q, the 
products of sums of independent excitations exponentiate: 

= expS (6) 

For all couplings < x < x c , and for general SU (N) and d, the exponential form (6) 
holds, where the function S is a single sum over the sites s of the lattice of all possible 
linearly independent "linked clusters" Cj a of SU (N) Wilson loops: 

s i 

By "linked" is meant products of Wilson loops at fixed relative separation and orientation. 
The "amplitudes" a% are then the coefficients of these clusters, summed to all orders 
in strong coupling perturbation theory. The representation (6), (7) holds, and is exact, 
provided the true vacuum is not orthogonal to the strong coupling vacuum ^>o(x = 0) 
i.e. up to the lowest x at which there is a phase transition. Thus, for any N and d, the 
Schrodinger equation for the vacuum energy is equivalent to 

E (E a (L) 2 S + E a (L)S ■ E a {L)S) -x^Re[/ p = ^» (8) 

L p s 

The wavefunctions of excited states (glueballs) with zero momentum can be written 

f 9 = Gex P S (9) 



where 

G ' = E( 6 o + E^ C 'm) (10) 

s i 

i.e. G is a similar sumP] over the lattice of linked clusters of SU(N) Wilson loops. The 
corresponding Schrodinger equation is then 

Y,(E a (LfG + 2E a (L)S ■ E a {L)G) = (w g - w )G (11) 

L 

The effect of the electric field operators on a given link L in (8), (11) is to recombine 
the clusters in S and G containing U(L) or U^{L) according to the relation for SU(N) 
generators T®pT® s = \8 a $8p^ — jj^8 a p8^s. Some examples for general SU(N) are 



where the dots indicate the ends of the link L. 

The exponential form (6) for the vacuum wavefunction is in fact a very general ex- 
pression for the fully interacting ground state of a many-body system, with S a sum of 
linked clusters each of which describes linked excitations of a bare or unperturbed state 
$?o(x = 0). The exponentiation then gives the correct statistical weighting for unlinked ex- 
citations consisting of products of linked components. Furthermore, it guarantees that the 
Schrodinger equation for the vacuum and for excited states is a single sum of linked terms 
only, so that there are then no size extensivity problems with approximation schemes. The 
expression (6) is the starting point for the many-body theory coupled cluster method^. 
This is an intrinsically non-perturbative approach, involving the amplitudes Oj,6j directly 
rather than their perturbative expansions. 

In order to apply the coupled cluster method to lattice SU(N) gauge theories, the 

fundamental approximation in general consists of two steps: 

3 Excited states of a given spin and parity are constructed through the appropriate dependence of the 

bi on the clusters' rotation/reflection orientations, left implicit in (10). 

4 For an introductory review of this method, which has been widely used in nuclear and condensed 

matter physics and quantum chemistry, see ref. 0. Application to lattice gauge theory was suggested by 

Greensite W, but this suggestion was not followed up. 
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I) A truncation of the sums of clusters in 5 and G according to the clusters' geometrical 
size. All clusters then formed in the Schrodinger equations with the truncated 5 and G 
which are larger than the given cut-off are simply discarded. 

II) A Taylor expansion of 5 and G in the remaining clusters, rearranged, or "shifted" , to 
describe fluctuations of the individual Wilson loops about their vacuum expectation values. 
These expectation values are calculated self-consistently from the vacuum wavefunction. 

The first step was originally suggested by Greensite but in order for the method 
to provide a tractable, valid and efficient approximation scheme beyond the perturbative 
range of couplings, the second step is in general essential. The physical idea is that 
amplitudes in 5 involving fluctuations on scales larger than the vacuum correlation length 
£, and in G larger than the size p of the glueball, are expected to be greatly reduced 
relative to those for fluctuations smaller than £ and p, and to decrease very rapidly with 
increasing size. For fixed lattice spacing, the results should therefore converge rapidly to 
the exact result as the dimensions of the cut-off cluster size are increased beyond £ and 
p. As the lattice spacing a — * 0, £/a and p/a diverge and the size in lattice units of 
clusters required to describe the theory must diverge correspondingly. The success of the 
method therefore depends on the results obtained at successsively increasing values of x 
(i.e. decreasing a) showing a sufficiently rapid approach to the expected scaling behaviour 
that the continuum limit can be extrapolated. 

We consider first the case of SU{2) with a one-plaquette cut-off, for which step I is 
exactly soluble. For 577(2), the identity TrATrB = TrAB + TrAB" 1 for 5*7(2) matrices 
A, B gives linear relations between clusters e.g. 



+ 



0= 2 +□ 



Thus, the one-plaquette clusters may be written as "powers" of the one-plaquette single 
loop Ci iP . Then, with the one-plaquette cut-off 5, writing the infinite sum of one-plaquette 
clusters on each p as a function S p (Ci tP ) 

S = J2S P (C liP ) (12) 
v 
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Schrodinger's equation for the vacuum becomes the differential equation 



where terms involving clusters spanning adjacent plaquettes p,p' 



have been discarded. The solution isQ 



dSp dSp, Eac ^ ^ (M) 



S p = In ^(1 - \Clp) Xl \e 2 {-\™s-\-\C x , v ) j j (15) 

where se 2 is a Mathieu function [pLOfl . The associated vacuum energy per site wo for <i = 3 
is plotted in fig. 1, together with Pade approximants obtained from the strong coupling 
series to 0(x 14 ) |llj with which it shows very remarkable agreement at all x. Furthermore, 
using the Feynman-Hellman theorem, the vacuum expectation value (C\) = —dwo/dx, has 
the correct weak coupling behaviour 

(d) = 2- 0{x- 1 / 2 ) x » 1 (16) 

The first excited state involves the Mathieu function se^. The mass gap Aro 9 = w g — wo 
for d = 3 is shown plotted against x -1 / 2 in fig. 2, together with Pade approximants from 
0(x 10 ) series [11]. In d = 3, the SU(2) theory is superrenormalizable and so the true 



(reduced) mass gap Aro 9 is expected to go to a constant at the continuum limit x — > oo. 

At weak coupling, the Mathieu mass gap diverges as Atu ff = 4s 1 / 2 — f + 0{x~ 1 / 2 ). This 

x 1//2 behaviour is as expected from the uncertainty principle with a wavefunction involving 

excitations spanning only a finite number of lattice spacings. 

In general, increasing the cut-off size beyond one plaquette or going beyond SU(2), 

the resulting partial differential equation for S p is intractable. We therefore introduce the 
5 Eqn. (15), which is the exact solution for an SU(2) "one plaquette universe", has been derived by 
many authors e.g. refs. ||, j?]], and also as the d — 3 variational solution il is taken to be a product 
of functions of one plaquette variables S]. However, the fact that it reproduces so well the results of ref. 
|lH for d = 3 has not been noticed previously. For d — 4 (for which a variational solution is not available) 
eqn. (15) still results in a vacuum energy which compares well with that obtained from Pade approximants 
derived from strong coupling series pi. 
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second step of the approximation scheme and rearrange the sums of clusters in the given 
cut-off S and G in terms of "shifted" linked clusters C[ s 

s = EE ffl ^ ( 17 ) 

s i 

G = EO'o + E^O (18) 

s i 

These shifted linked clusters consist of products of Wilson loops each minus their vacuum 
expectation values, still at fixed relative separation and orientation e.g. for general SU(N) 

(D)' with = □ - (□>, □' = □ - (□) 

To illustrate the method, we consider again the one-plaquette cut-off and include in S 
and G shifted one-plaquette clusters up to the n'th power of the shifted one-plaquette loop 
C' lp = C\ tP — (C\). Substituting into (8) and (11), the resulting clusters spanning two 
plaquettes are discarded (step I) together with all "higher order" fluctuations (C' l p ) m , 
m > n, (step II) and coefficients of the constant term and each shifted cluster (C{ p )*, 
i = 1 to n, are equated. For the vacuum, this gives a set of n non-linear equations for the 
amplitudes a[ and an expression for wo, all in terms of x and (C\). For the glueball, the 
procedure gives an n x n matrix equation for the b[ with eigenvalues Aw g , together with 
an equation for b' , all in terms of the a- and (Ci). These coupled cluster equations are 
then solved using numerical library routines, with (Ci) calculated self-consistently 



The results for wq and Aw 9 obtained |12[ with the one plaquette cut-off and the 



expansions of S and G truncated at n = 1, 2, 3 . . . converge rapidly to the Mathieu results 
at all x. Each successive power of C[ p generates the next term in the series expansions 
in x for x ^ 1 and x -1 / 2 for x ^> 1 of the Mathieu eigenvectors and eigenvalues. In 
particular, even at the very simplest n = 1 approximation the method gives the correct 
weak coupling leading behaviour (16) for the vacuum expectation value (C\). 

Proceeding to larger cut-offs with SU(2) in d = 3, we have included in S and G all 
the shifted clusters which, unshifted, occur in 2nd, 3rd and 4th order strong coupling 
perturbation theory (which is used solely as a guide to which terms to include - the 
resulting coupled cluster equations are all solved non-perturbatively). The complete set 



x order 1 order 2 order 3 order 4 Pade |11] 



0.5 


-0.106298 


-0.0843647 


-0.0827613 


-0.0822829 


-0.082156 


1.0 


-0.386012 


-0.332963 


-0.322072 


-0.319032 


-0.31587 


1.5 


-0.779896 


-0.707363 


-0.685324 


-0.681253 


-0.67115 


2.0 


-1.25000 


-1.16553 


-1.13371 


-1.13132 


-1.1162 (2) 


2.5 


-1.77435 


-1.68175 


-1.64112 


-1.64197 


-1.626 


3.0 


-2.33949 


-2.24077 


-2.19185 


-2.19667 


-2.183 (3) 


3.5 


-2.93665 


-2.83296 


-2.77600 


-2.78509 


-2.77 


4.0 


-3.55974 


-3.45180 


-3.38692 


-3.40037 


-3.40 (2) 


4.5 


-4.20439 


-4.09265 


-4.01992 


-4.03771 




5.0 


-4.86730 


-4.75207 


-4.67154 


-4.69358 




10.0 


-12.1254 


-11.9822 


-11.8255 


-11.8828 




20.0 


-28.2197 


-28.0321 


-27.7362 


-27.8358 





Table 1. Vacuum energies Eo for the first four orders of approximation in the shifted coupled cluster 
method, together with the results of Pade approximants derived from the 0(x 14 ) strong coupling series 
taken from fnl] (estimated errors in the last figure are given in brackets). 



of 69 linearly independent clusters occuring at 4th order are shown in the appendix. The 
loop expectation values are calculated self-consistently from the wavefunction using the 
Feynman-Hellman theorem by making the change 

W^W + Y, E e iCi,s (19) 
s loops eS,G 

so that the a- become functions of the ej and x. Then 

dw 



(Ci. 



(20) 

fe=o> 



de 

These are then simple quantities to calculate numerically, avoiding any group integrals 



The results, obtained using a computer to carry out all the calculations [12], are 
given in table 1 for wq and shown plotted in fig. 3 for Aw g . Assuming that the Pade 
approximants are a good guide to the exact results, the method gives excellent results 
for the vacuum energy wq at strong and intermediate couplings and the correct form 
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wq = —2x + 0(x 1 ^ 2 ) at weak coupling. All the Wilson loop vacuum expectation values 
have the correct weak coupling form, as in (16). For the mass gap Aw g , reasonable results 
are obtained at intermediate couplings, before the inevitable x 1//2 scaling sets in at weak 
coupling. In addition to the eigenvalues, the method also provides the wavefunctions 
and the cluster vacuum expectation values at all couplings. The amplitudes Gtj, bi of 
the most dominant of the 69 shifted clusters occuring in the wavefunctions at 4th order 
are shown in figs. 4 and 5 respectively. At weak coupling all the amplitudes behave as 
a!. = a- (1/2 V/ 2 + C>(const.), b\ = ^ (0) + C(x" 1 / 2 ), with, for the 69 clusters occuring at 4th 
order, the a'^ 1 ^ 2 ^ and b^ ' ranging over six orders of magnitude. For the vacuum, it is found 
that, at weak coupling, simple shifted clusters spanning few plaquettes and/or containing 
few Wilson loops dominate over those for complicated multi-loop clusters. This can be 
understood as a combinatorical effect in the coupled cluster equations. For the lowest 
(scalar) glueball, shifted clusters consisting of single Wilson loops dominate at intermediate 
and weak coupling. With the glueball calculation involving the diagonalization of a (large) 
matrix, this effect is harder to understand. No evidence is found of a phase transition at 
finite x. 



Returning to general SU(N), it can be shown [13] that the shifting procedure guar- 
antees that, for any N and d, the approximation scheme gives the correct weak coupling 
vacuum expectation value N + 0(x~ 1 ' 2 ) for a Wilson loop in S. Thus, at x — > oo, the 
effect of the electric field operators coupling together e.g. a pair of shifted single loops 
becomes 

E a (L)D!-E a (L)j-^i =I^g'-D'- yj Q' + const. term 0(x- 1/2 ) 

It is straightforward to show that, in general, in the weak coupling, large N limit, only 
shifted single loops[] survive in S and G. Why this appears to be a good approximation 
at N = 2 remains unclear. 



6 For SU(N), a single Wilson loop can "wrap around" itself up to N — 1 times before it becomes linearly 
dependent on clusters of simpler loops. 
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Conclusions 

In general, the application of the coupled cluster method to lattice gauge theories 
depends crucially on the rearrangement, or shifting, of the clusters in S and G to describe 
gauge invariant fluctuations. At x <C 1 the method matches on to perturbation theory^, 
while at x S> 1, it is always found that small, simple shifted clusters dominate in the 
wavefunctions over large, complicated shifted clusters due to combinatorical effects in the 
Schrodinger equations. Furthermore, the interpretation of the discarded terms in the 
Schrodinger equations as higher order fluctuations enables otherwise formidable problems 
involving linear dependences among clusters to be largely circumvented. 

In summary, the shifted coupled cluster method provides a non-perturbative semi- 
analytic approach to lattice gauge theories, involving a direct physical approximation 
scheme, and gives non-perturbative information on both eigenvalues and eigenfunctions 
of the Schrodinger equations. The results presented here show (table 1) that the method 
converges very rapidly for the vacuum energy (assuming that the Pade results are a good 
guide) and should be capable of providing excellent approximations to the vacuum wave- 
function. Further work is needed to establish how well it works for excited states; through 
a more judicious choice of clusters in S and G, it should be possible to improve signifi- 
cantly upon the results presented here for the mass gap at intermediate couplings so that 
it can be reliably extrapolated to the weak coupling limit. The method therefore promises 
to provide much valuable analytic information on the behaviour of lattice gauge theories. 



Fuller accounts of this work will appear elsewhere [13]. 
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Appendix: Fourth Order (Unshifted) Clusters 



Cl,s — 


□ 


<^2,s = 


I II 1 

□□ 


c = 






O 


C5,s = 


□ 


C§ yS = 


□ 

1 1 


CV,s = 


on 


C8,s = 


!□ 1 


C9,s = 


□□□ 


ClO,s = 




Wl,s — 


jn 


Cl2,s = 




Cl3,s = 





Cl4,s = 


□ 

□□□ 


C15 s = 


□ 



Cl6,s = 

C\7,s = 

Cl8,s = 

Cl9,s = 

C20,s = 

C21, s = 

C22,s = 

C23,s = 

C24,s = 

C25,s = 

C26,s = 

C27,s = 

C28,s = 

C29,s = 

C30,s = 

C31,s = 

C32,s = 

C33 . = 



□ 

on 

□□ 
□□ 



□ 
□0 



□o 



□ 



□ 



□ 

□ 
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□□□ 

□ 



□ 



□ 



C34,s = 
C35,s = 
C36,s = 
C37,s = 
C38,s = 
C39,s = 

Cao, s = 

C41,s = 

C42,s = 

C43,s = 

C44,s = 

C45,s = 

Ca&,s = 

Cw,s = 

C4S,s = 

C49,s = 

C5Q,s = 

C51 « = 



□Dn 
da 



□ 



□□D 



ODD 



□on 



□ 



□ 



□ 



□ 
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Cg2,s = 
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Cq4,s = 

C&5,s = 
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□ □ 
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Figure Captions 

Fig. 1. Graph of the vacuum energy wq vs. x for the exact Mathieu solution of the 
one-plaquette-cutoff coupled cluster equation (13), together with i) the Pade approximants 
derived from the 0(x ) strong coupling perturbation theory series taken from [11], and 
ii) the results from 0(x 2 ) and 0(x 4 ) strong coupling perturbation theory, which show that 
the Mathieu result is non-trivial. 

Fig. 2. Graph of the mass gap Au> 9 vs. x~ 1 / 2 for the exact Mathieu solution of 
the one-plaquette-cutoff coupled cluster equation, together with i) the Pade approximants 



derived from the 0(x ) strong coupling perturbation theory series taken from [11], and 



ii) the estimated bounds, shown dotted, for their extrapolation to the continuum limit 



x — > oo, also from pj |. 



Fig. 3. Graph of the mass gap Au> 9 vs. x -1 / 2 for the shifted coupled cluster wave- 
functions fy g = GexpS with S, G including the sets of clusters which, unshifted, occur 
at the first four orders of strong coupling perturbation theory. Also shown are the Pade 
approximants and continuum extrapolations derived from the 0(x 10 ) strong coupling per- 
turbation theory series taken from [11]. 

Fig. 4. Graph of the absolute values |a^/a^ | vs. x for the amplitudes a[ of the most 
dominant of the 69 shifted clusters occuring in the vacuum wavefunction = expS" at 
the 4th order approximation, where is the amplitude of the shifted single plaquette loop 
C[ . The cluster labels i are shown on the right side of the figure. 

Fig. 5. Graph of the ampitudes 6^ vs. x for the most dominant of the 69 shifted clusters 
occuring in the scalar glueball wavefunction ty g = GexpS at the 4th order approximation. 
The cluster labels i are shown on the right side of the figure. 
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